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Wave  Propagation  in  Laterally  Varying  Media:  A  Modal  Expansion  Method 

by 

Charles  B.  Archambeau 


Zonal  Partitioning  and  Green’s  Function  Representations 

Consider  a  two  dimensionally  varying  elastic-anelastic  medium,  as  indicated  in  Figure  1 . 
In  each  zone  Vj,,  a=  1,2,  ••  •  M,  the  medium  varies  in  the  vertical  direcnon  (z),  but  is  uniform  in 
the  horizontal  direction  (y  or  p).  The  supposition  is  that  the  laterally  varying  medium  can  be 
approximated  by  a  serit  s  of  step  variations  in  material  propeities  in  the  same  way  as  is  done  in 
the  vertical  direction. 

In  Vg  wa  have  for  he  frequency  domain  displacement  field  at  any  point  r  within  Ve,: 


(«)u,(r,  (D)  i 


_L 

47t 


J  tj(ro)(«’G/(r,  ro;  (o)  -  Uj(ro)(“)g/(r,  ro;  w)! 

£o  +  £a-l'-  ■’ 


(1) 


where  and  (“^g/  are  the  zonal  displacement  and  traction  Greens’  functions  appropriate  for 
the  zone  or  region  Va  *  T'he  vertical  boundary  surfaces  of  are  la  and  Ea-b  ^  indicated  in 
Figure  1.  Here  we  assume  no  sources  iiirid-"  and  that  the  Green’s  functions  satisfy  all  inter¬ 
nal  boundary  conditions  on  all  horizontal  layers  in  Va.  (In  this  case  there  are  no  surface 
integrals  over  intei  lal  boundaries  in  (1)).  Green’s  functions  in  Va  can  be  written  in  lerms  of  the 
eigenvalues  ka  and  eigenfunctions  (“V  for  this  zone  as  ** : 


f  • 

Summdion  over  repeeted  coordmete  indicei  is  used  throughout.  Coordinate  ind'ces  wiU  appear  as  lower  case  latin  sub- 
smpta  and  supersaipts.  The  summation  convention  does  nui  apply  lo  any  indices  appearing  in  parenthesis. 

Throughout  this  development  the  "sum "  over  the  eigenvalues  ka  will  be  written  as  a  discrete  summation  but  it  should 
be  understood  that  ui  an  unbounded  medium,  sucli  as  a  layered  half  space,  part  of  the  wave  number  spectrum  will  be  con¬ 
tinuous  In  this  case  the  "sum "  over  ka  must  be  interpreted  as  a  generalized  summation  involving  a  regular  sum  over  the 
disci,  le  part  of  the  spectrum  plus  an  integration  c  er  the  continuous  pan  of  the  wave  number  spectrum 
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5)a-i  2a  2a+i 

I  i  i 


Figure  !•  Zonal  partitioning  of  a  vertically  and  laterally 
varying  medium  into  subregions  of  uniform 
haizontal  layering. 
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(2) 


To;  co)  =  4n  2 

m.kfl 


kg)  ^°Vi(r.  kg) 

Na(^a  > 


where  is  the  complex  conjugate  of  ^“Vi  and  No  is  a  normalization  constant  which  may  be  a 
function  of  frequency  w  and  the  wave  number  ka.  Since  the  are  eigenfunctions  for  the 
region  Vc  this  Green’s  function  satisfies  all  boundary  conditions  along  the  horizontal  boun¬ 
daries  in  Va.  (For  details  see  Harvey,  1983.) 

Further,  since: 

To;  ©)  =  ni[o) 

where  n|^°)  is  the  surface  normal  to  Ea  and  Ea_i  and  x/°)  are  source  coordinate  variables,  then 


a(«)G^(r,  To;  (£>) 
ax/o) 


(«)g/(r;ro;  tD)=:4n£ 

m,k,i 


(ro.  ka)Wn/.(r.kc) 
Ro^lcoT'to) 


(5) 


Here  m  is  the  angular  index  for  cylindncal  coordinates,  ka  the  horizontal  wave  number 
corresponding  to  the  modes  in  Va  and  where: 


(«)H5(ro,ka)  =  n;[o)a{})n 


Bxf 


^“Vn(ro.  ka) 


(4) 


Because  of  the  horizontal  layering  in  Vq,  the  eigenfunctions  and  are  defined  sec- 
tionally,  that  is: 


|(a)\j/^(s)(z) 


Zs_i  ^  Z  ^  Zs 


with  (s)  the  honzontal  layer  index  in  Vq. 

For  the  horizontally  layered  region  Va  we  have  that: 


(«)G/  =  +  («JG/I 

(«)g/=(«H'Rg/  +  'fg/ 


(5) 
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Here  arid  G/  are  the  Rayleigh  and  Love  type  Green’s  displacement  funct.ons  (with  simi¬ 
lar  names  for  the  associated  Green’s  tractions)  and  wheie 


faV'i/-  -  .  V*  Lka)^®^V|^(ro,  I  ka) 

(t)0/(r.r.,co)  =  4;._L  ^  o) - 


with  Rka  and  tka  representing  the  Rayleigh  and  Love  type  mode  eig  nvalues.  Likewis 


In  cylindrical  coordinates  (p,  (ji,  z),  the  eigenfunctions  are  (see.  for  example,  Harvey,  1981;: 

N 

<^V(r.  ska)  =  C);n(z;  Rkg)  P(Rkap. 

+  («)Em(z;Rka)B(RkaP.(l))i 
(“R)'P(r,  Rka)  =  («)R„,(z;  Rka)  Pm(RkaP.<|)) 

^**^Sm(z;  Rka)  B[n(Rka  p,  <})) 


^"^(r.Lka)  =  ^‘*^Fn,(z;  tka)  Cjtkap,  <|)) 
^(r,  Lk,)  =  («)7’n,(z;  Lka)Cn,(LkaP.(t)) ' 

Here  P,  B  and  C  are  the  vector  cylindrical  harmonics  defined  as: 


Pm(kp,(l))  =  ejJn,(kp)e“"* 

B„(kp,«.  6p^+e,  ^  J„(kp)e»'* 


Cn,(kp,(l))=  Cpl-j^  Ukp)e'''’*! 


where 
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Jm(kp)  =  H;i,»)(kp)+H^2)(kp) 


with  the  cylindrica!  Bessei  function  and  and  the  cylindrical  Hanke!  functions. 
These  vector  functions  are  clearly  such  that  Pm’Bm  *  Pm  Cm  *  Bm  Cm  =  0  and  also  have  the 
usual  functional  orthogonality,  (e.g.  Stratton  1941.  Morse  and  Feshbach.  1953).  Here  ii,  ip  and 
are  the  unit  vectors  in  cylindrical  coordinates.  The  various  "stress-displacement"  functions 

in  (8)  are  the  same  as  those  usually  appearing  in  the  ordinary 
developments  for  a  laterally  homogeneous  layered  half  space  -  such  as  described  in  Harkrider 
(1964);  Ben  Menahem  and  Singh  (1972),  or  Harvey  (1981). 

Similar  representations  for  the  eigenfunctions  can  be  given  in  cartesian  and  spherical  coor¬ 
dinates.  (In  the  latter  case  the  eigenfunctions  kH'  and  lHi  are  usually  termed  spheroidal  and  tor¬ 
sional;  and  P,  B  and  C  become  vector  spherical  harmonics).  The  choice  of  cylindrical  coordi¬ 
nates  implies  rotational  symmetry,  that  is  that  the  medium  is  partitioned  into  zones  Va  which  are 
cylindrical  shells,  with  "igure  1  depicting  a  cross  section  at  lixed  (j>.  If  cartesian  coordinates  are 
used,  then  Figure  1  represents  a  cross  section  at  constant  y,  with  properties  constant  in  the  ±y 
directions.  In  the  development  that  immediately  follows  cylindrical  coordinates  will  be  used; 
however  the  cartesian  and  spherical  representations  are  also  appropnate  and  the  development 
and  results  are  analogous  to  tliose  for  the  cylindrical  choice. 

"Forvrard"  and  "Backward  Propagating"  Mode  Expansions 

In  addition  to  the  eigenfunction  expansions  of  the  Green’s  functions  in  Vq,  we  can  also 
expand  the  displacements  and  tractions,  appeanng  in  (1)  in  terms  of  eigenfunctions  in  Vg.  In 
particular,  ,(ro)  and  ^“Hj(ro)  may  be  expanded  in  terms  of  "forward"  and  "backward"  pro¬ 
pagating  modes  as: 
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<®^Uj(ro.  o) « <®>u/”(ro.  co)  +  <®^u/^Kro.  co) 

(*Hj(ro.  w)  a  »(ro,  (o)  +  <o) .  ^  ^  ^  ^ 

where  the  superscripts  (1)  and  (2)  denote  modes  pr(^agatmg  in  the  positive  and  negative  radial 
(p)  directions.  Specifically, 

r  1 

w'Jii*- 

where 

<®y»»>(roJ.k;)  a  [<«)Dm^Zo:  K)  P^P>  +  <“>E„-(Zo;  lOB^  +  <®)F„.(Zo;  K)  C^p)]  e'"'^ 

(13) 

(®W(ro.  ki) .  [w/?„,.(ao;  ki)  +  <«5S„.(Zo:  k^)  +  r„.(z<,:  k;)  CMi*']  e‘«"> 

with 

» 

•  L  J  J 

J 

The  coefficients  are  to  be  determined  from  boundary  conditions  at  £«  and  la.i, 

these  conditions  bring  the  continuity  of  displacement  and  traction  on  these  surfaces.  On  the 
other  hand,  of  course,  all  the  functions  and  are 

known  functions  of  the  coordinate  variables  and  the  intnnsic  matenal  properties  of  tl  2  internal 
horizontal  layers,  since  they  are  provided  by  the  usual  one-dimensional  propagator  approach  in  a 
layered  half  space  (e.g..  Harvey,  1981).  The  explicit  forms  of  the  functions  are  included  in  the 
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Appendix  1. 


Given  that  0}  and  g/  in  (1)  can  be  split  into  Rayleigh  and  Love  type  Green’s  functions,  as 
defined  in  (5)-(9).  then  it  follows  that  <®^Uj  can  also  be  split  into  modal  sums  involving  only 
and  Therefore: 


<“>Uj  a 


(15) 


where 


P"l  m*.  Dk^ 


E<fu/^Hro.w)«  E 

p"l 


(2L 


(16) 


with 


<«)D^.(Zo ;  «lCa)P#  +  <»)E„,.(Zo :  Hk„)  e^^'* 


WP>(ro.Rk«)i 
<?>'  ‘Hto,  ika) « <“^F„.(Zo ;  iJca)  C^|))  e“"'* 


(17) 


A  similar  decompi  >sition  applies  to  the  u  action 

It  is  importart  to  note  that  the  eigenfunctions  used  to  e  xpand  the  Green’s  functions  in  equa* 
tions  (2)  -  (7)  are  appropriate  for  the  horizontally  layered  zone  m  Va  and  are  themselves  normal¬ 
ized  such  that: 


<(«^V/P)(ka).  WPKk;)>  a  I  <“R>v|//P)(kar)(®K>;ii^P)(kir)dV  =  6(k„  -ki)  6ff' 
<(«LV/PHka).  =  J  f ^?V/PHk;r)dV  =  6(ka  -kX 


(19) 


where  denotes  the  complex  conjugate  of  and  the  right  hand  side  involves  the  usual 


delta  functions.  Therefore  the  nonnalization  factors  appearing  in  the  Green’s  function  expan¬ 
sions  are  free  parameters  that  may  be  chosen  so  as  to  appropnately  normalize  the  zonal  Green’s 
functions  in  Vn,  a  =  1,2,  •  •  •  M. 


To  obtain  the  ^propriate  normalization  factors  for  and  and,  in  addition,  to 
express  these  Green’s  functions  in  forms  that  are  convenient  for  use  with  the  expanded  form  for 
(®)Uj  in  (15)-(16),  it  is  useful  to  adopt  an  expansion  form  for  the  Green’s  functions  that  is  similar 
to  that  for  (“^uj  in  (15).  That  is,  using  both  and  in  the  expansion  for  (“>G/,  we 
express  the  Green’s  functions  as: 

N 

(«r>G/  (r.  To;  <o)  = 

G/  (r,  To;  (0)  =  (f  G,y)  +  )( 


where: 


ralQ (oI  -  V  Rk  .) (“K)V,(P)(r,  Rka) 

_  RNnk^T^i^)  ! 

(«)G(P)-4ny  LkJ  i 

nTt  LN^“>(ka,a))  I 


and  similarly  for  and  the  Green’s  tracnons. 


Orthogonality  and  Normalization  Relations  for  Zonal  Eigenfunctions 

We  can  use  (15)-(16)  in  (1)  and  also  substitute  (20)-(21)  into  this  representation  integral. 
Since  the  representation  given  by  (1)  should  be  of  the  form  of  the  expansion  in  115),  we  should 
obtain  by  proper  choice  of  the  normalizanon  factors,  and  exactly  the  expansion 
given  in  (15)  in  terms  of  forward  and  backward  propagating  modes.  In  particular,  from  ( 1 )  we 
have: 

l“)uj  (r,  (o)  =  l*^Uj  (r,  co)  +  (^^Uj  (r,  co)  (22) 
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with 


(•^UiCr,  CO)  = 


ih 


•RUj 


1 


(>,{r,co)=  f  Ltj  (fg; 


dflo  ;«’eVa 
dflo  ;reVa 


(23) 


Introducing  the  exphcit  eigenfunction  expansions  from  (15)-(16)  and  (20)-(21),  we  get: 


f<^Uj(r,  (o)  =  co) 

p“i 


(24) 


m'.Rtci  m.Rka  rN^  '  (m  .  W) 


<'«PKRk;).  W*>(Rk„)>„,„.i 


-  <(“R'H'/P)(Rk;).  WKRka)>a,«-l 


-  <  (Kk;),  (Rka)>„.  (r,  k  J 


(25) 


Here  terms  of  the  form: 


^M^j(ka)'  Xj(ka)^a,  o-l  *  ^'l^j(ka)'>  Xj(ka)^a  "f"  ^Vj(ka)i  Xj(ka)^ a- 


are  introduced,  where  the  inner  product  is  defined  over  the  surface  La  (or  la-i)  as: 


<Vj(k;),  Xj(ka)>  s  ‘^«^X)('’o’  ka)dao 


with  summation  over  the  repeated  coordinate  index  (j)  implied.  An  exactly  analogous  result 
holds  for  lUjI  with  the  suffix  "R"  replaced  by  "L"  in  (24)  and  (25), 

Comparing  (25)  with  the  equivalent  expressions  in  (15)  -  (16),  it  is  clear  that  the  mner  pro¬ 
ducts  appeanng  in  (25)  must  reduce  to  delta  functions  over  the  angular  index  m  and  the  mode 
eigenvalues  ka-  In  particular,  the  following  orthogonality  conditions  apply*: 

'Where  it  is  obvious  Croin  oontext,  the  R  and  L  identifying  subscripts  on  the  wave  numbers  Rkg  and  Lkg  will  be 
suppressed  in  order  to  reduce  clutter  in  the  e<)uations. 
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<(“^'F/P)(k;)(«^vj//P)(ka)>p  -  <(«^v/pHk;),  Wp)(ka)>| 


=ii 


('^'F(p)(k;ro)  •  <‘S)viJ(P)(karo)  -  <'^VP’(k;ro)  •  (‘^VP)(karo) 


d^o 


(26) 


=  np  [Pip>(k' pp)  •  P^fKkaPp)  +  B^Hk^p)  •  B^)(kaPp) 


with  P  =  a,  a  -  1  and  p  =  1,2  and  where  np  =  2npp.  (1  !ere  pp  is  the  constant  value  of  the  radial 
coordinate  on  the  surface  £p.)  In  addition, 


<(««>'Pj(P)(k;).  (“^v/‘J)(ka)>p-<(“K>WPHk;).  (“R>'P/‘»(ka)>p  = 

I  [(“K>'P/P)(k;ro)  •  (“^?‘i)(k«ro)  -  WPHkiro)  •  (“R>V‘i){karo)]  d^o  =  0 


(27) 


for  P  =  a,  a  -  1  and  p  ^  q.  Formally  identical  relations  hold  for  the  eigenfunctions  (‘[)y  P)  and 
are  obtained  by  replacing  the  suffix  "R"  by  "L"  in  (26)  and  (27).  Here  we  observe  th  it  the  for¬ 
ward  and  backward  propagating  modes  are  completely  orthogonal  sets.  These  conditions  are 
equivalent  to  those  obtained  by  Herrera  (1964)  and  McGarr  and  Alsop  (1967)  and  were  used  by 
Kennett  (1983)  in  his  development  of  a  formalism  for  wave  propagation  in  laterally  va^'ing 
media.  In  more  explicit  form,  equations  (26)  and  (27)  are  equivalent  to; 


|[(«)R„(Zo ;  k;)  («)I^(zo  ;  k«)  -  («)D„,(Zo ;  k^)  ;  k„)|dZo  = 

If “^S„,(Zo ;  k;)  («)E„(zo ;  k^)  -  <«>£ Jzo ;  k^)  Zo ;  kjj  dz^  = 


where  the  subscript  "R"  on  the  P-SV  wave  number  has  also  b^  en  suppressed  n  these  expres¬ 
sions.  For  the  SH  modes  the  analogous  orthogonality  relanon  is  easily  seen  to  be 


|[(«)T^(Zo  ;  k;)(«)Fn,(Zo  ;  k„)  -  (z^ ;  k^)(«)T^(Zo  ;  k«)  i  dz^  =  6.'  ;■ 
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where  the  wave  numbers  and  and  now  refer  to  the  SH  wave  number  set  Lko,.  The  "vertical 
eigenfunctions"  in  Va  are  those  defined  in  (13)  and  are  simple  expontials  in  Zq.  (See  Harvey, 
1981.)  Here  also  we  consider  the  kg  to  be  discrete  infinite  sets,  so  that  orthogonality  is  expressed 

by  the  Kronecker  delta  6^. 

Using  these  orthogonality  relations  in  (25)  gives: 


(«R)u/»(r,a))=  £  <faff(ka)Z 

m',  k<  m,  ko 


RN(“>(ka,  CO) 


nalPiii^Hk;pa)  •  Pi^'KkaPa)  +  B,(,'')(k;pa)  •  B^i‘>(kapa) 


+  na-iiP,!,'  (k;pa-i)  •  Pi(i‘HkaPa-l)  +  B,^^^(k;pa.i)  '  B^”(kaPa-l) 


ka) 


So 


(«R)u/n(r.co)=  £  («K>ai,‘^(k«)  («R'v/‘Hr.ka) 

m,  ku 


provided  we  take: 


(kaPa)  ■■  P,i,‘KkaPa)  +  B^»  (kaP„)  '  B^‘)(kaPa)| 


+  na-liP,^'Kk«pa-l)  ■  P^i‘HkaPa-l)  +  B^‘)(kaPa-l)B^,‘KkaPa-l) 


(28) 


Similarly, 


Ru/2)(r,co)=  £  («K>«m\k«)W>(r,ka) 


m.ka 


provided 


= 


na  P^2)(k^pj  .  p^2)(k«p„)  +  BA?HkaPa)  '  B,^2)(kaPa)  S' 


(29) 
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+  n«.i|p^2)(kap«_,)  •  P^2)(k^p^_^)  +  B(?)(k„Pa.,)  ■  B«)(k«p„_,)[ 

Tlie  results  for  are  analogous  and  the  normalization  factors  are: 

=  [naC^l)(k„Pa)-C^‘>(kaPa)  +  Ha-lC  (,'HkaPa-l)^'Hk„p«.,)]  (30) 

lN2(«^  =  naC^2)(l^^^)  .  +  na_|Ci2)(l^^p^_l)^2)(k^p^_l  j  (31) 

Thus,  the  form  of  the  displacement  field  in  any  one  of  the  zones  Vq  is  given  by 

(«Wr,t0)=  £  [(“R>ai‘\Rka)W‘H^Rka)+(“R>Om^(Rka)(«M^  <ka) 
m.Rka'- 

+  E  i  r  '  V„  (32) 

m.LlCa'-  ■‘ 

which  is  (merely)  a  sum  of  P-SV  modes  propaganng  in  the  forward  and  backward  honzontal 
directions,  plus  a  similar  sum  of  SH  modes.  Further,  the  displacement  field  in  V^,  is  connected  to 
its  values  on  the  boundary  surfaces  Ig  and  la-i  hy  the  representations  in  (23),  with  the  Greens 
functions  given  by  the  eigenfunction  expansions  of  (20)  -  (21)  and  with  the  normalizations 
specified  by  (28)  -  (31).  Use  of  these  latter  representanons  provide  the  means  of  determining  the 
coefficients  and  in  (32),  and  thereby  an  explicit  expression  of  the  displacement 

field  in  Vo  in  terms  of  the  modes  of  this  horizontal!)  layered  region.  As  will  be  shown,  tlie 
coefficients  between  all  the  zones  Vq,  a  =  1,2,  ■  •  ■  ,  M,  are  linked  by  a  propagator  formalism. 

Zonal  Boundary  Conditions,  Projections  and  Lateral  Propagators 

Continuity  conditions  expressing  conseiwation  of  momentum,  mass  and  energy  apply 
throughout  the  medium,  however  complex  the  intrinsic  material  propeities.  In  particulai  such 
conditions  apply  along  the  control  surfaces  Ea  separating  the  zones  of  uniform  lateral  piopenies 
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in  Figure  1.  In  the  case  of  a  solid  medium,  with  welded  contacts  at  all  layer  boundaries,  the  con¬ 
tinuity  conditions  along  the  surface  are: 


'  ' 

t«>Uj 

(“Hj 

a 

«  . 

(33) 


where  the  subscript  a  on  the  matrix  brackets  is  used  to  indicate  evaluation  on  the  vertical  boun¬ 
dary  le,  between  the  zones  and  Vg+i 


The  displacements  and  tractions  in  (33)  can  be  expressed  in  terms  of  the  eigenfunction 
expansion  of  (32).  However,  smce  the  P-SV  and  SH  modes  are  decoupled  in  Va  and  1,  then 
(33)  can  also  be  expressed  by  the  decoupled  set  of  relations: 


LLt 

m  rIco  p»l 


m  Rka.i  P"l 


(«+‘)RA^>(Rka+i) 


(34a) 


LLt 

lICo  P“1 


m 


=E  E  t 

1^*1 


m 


<“^‘)LA^f)(tka+,) 


(“*lWj(P)(Lka+l) 


;j  =  3 


(34b) 


where  the  expansions  in  P-SV  and  SH  t  loves  have  been  substituted  for  Uj  and  tj  on  both  sides  of 

(33) .  A  similar  set  of  boundary  equations  apply  to  the  other  vertical  boundary  of  Vq,  on  the  sur¬ 
face  £a_  1.  in  Figure  1.  (In  this  case  the  matnces  are  evaluated  on  £«- 1  so  the  matrix  indices  in 

(34)  change  to  (a  -  1)  throughout,  while  on  the  right  side  of  (34)  all  the  eigenvalue  and  eigen¬ 
function  indices  change  from  a  +  1)  to  (a  -  1).) 

We  can  extract  expressions  for  individual  mode  coefficients  and  appropri¬ 

ate  to  the  zone  Vq.  in  terms  of  the  mode  coefficients  in  the  zone  Wa*  \  by  taking  integral  inner 
products  ("projecrons")  between  the  displacement  and  traction  eigenfunctions  on  both  sides  of 
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(34).  Then  we  can  use  the  (P-SV)  orthogonality  relations  in  (26)  -  (27),  along  with  comparable 
orthogonal  relations  for  SH  modes.  Specifically,  using  inner  product  bracket  notation  as  before 
in  equation  (25)  and  taking  the  inner  products  between  displacement  and  traction  eigenfunctions 


on  both  sides  of  (34),  we  have: 


m.ka  p«l 


<(«)V)//P)(ka)  ,  («)'P/»Hki"))>« 
<(«)'P/P)(ka)  ,  («)v|/j(*)(ki/»)>a 


=  E  t 

m, ka  p“l 


(«+')A(f)(k„.„) 


<(«+>  'V/PHka+i )  ,  '«)T/*'(k('’))>a 
<(“*l^'Pj(P)(ka+i)  ,  '«V/''(k^'’))>a 


(35) 


where  indices  R  or  L  have  been  suppressed  but  are  implied,  with  appropnate  use  depending  on 
whether  j  =  1, 2  or  j  *  3,  as  indicated  by  (34a)  and  (341)).  (That  is,  this  equation  applies  to  either 
(34a)  or  (34b)).  For  specificity,  one  uses  P-SV  eigenfunctions  and  eigenvalues  and  a  subscript 
"R"  when  considering  component  equations  with  j  »  1  2  and  uses  SH  eigenfunctions  and  eigen¬ 
values  with  subscript  "L"  when  considering  the  j  *  3  component  equation.)  Here  k^"'  denotes  the 
specific  n'**  eigenvalue  of  one  particular  mode  with  angular  index  m'. 

Now  we  can  subtract  the  upper  matnx  equation  in  (35)  from  the  lower  one  and  then  make 
use  of  the  orthogonality  relations  for  P-SV  modes  in  (26)  -  (27),  and  the  obvious  similar  pair  for 
the  SH  modes,  to  obtain; 


i  L  1  (“*‘)A^Kka+l)  <(«*‘V/PHka+l)  <“>'P/»>(k4'»)>a  - 

ka.l  P«1  I 

<(a+‘^'P/PHk„+i) .  ;  S  =  1.  2  (.36) 


where  we  have  equated  the  sums  over  m,  on  each  side  of  (35),  term  by  term.  Fhis  equation 
again  applies  to  either  P-SV  or  SH  modes;  however,  for  P-SV  modes  j  =  1,  2  and  or  SH  modes, 
then  j  =  3.  Therefore  in  (36)  the  implied  summation  over  the  coordinate  index  is  )ver  j  =  1  and 
2,  for  the  P-SV  case,  and  for  SH  modes  only  the  one  term,  for  wh  oh  j  =  3.  oc  urs.  The  fiee 
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index  (s)  denotes  the  forward  and  backward  horizontally  propagating  modes,  so  that  (36) 
expresses  a  relationship  for  both  mode  types.  The  factor  is  the  normalization  "constant" 
appropnate  for  the  different  mode  types.  These  factors  are  given  in  (28)  -  (29),  for  the  forward 
and  backward  propagating  P-SV  modes,  and  in  (30)  -  (31)  for  the  SH  modes. 

It  can  be  seen  from  (36)  that  a  particular  mode  in  V,,,  at  a  particular  eigenvalue  (or  wave 
number),  will  be  "excited"  by  all  the  forward  and  backward  propagating  modes  in  Va+i  in  the 
manner  descn  led  by  the  expression  on  the  right  side  in  (36).  Thus,  all  the  modes  in  Va+|.  at  all 
wave  number ,  will  contribute  to  the  excitation  of  any  one  mode  in  Va  (at  a  particular  wave 
number)  in  prt  portion  to  the  sum  of  tiie  mode  coefficients,  i®+'^A;^)(ka+i).  weighted  by  the  inner 
product  factors  given  by  the  bracket  teim  on  the  right  side  of  (36).  Thus  the  weight  factors  in 
(36)  will  be  called  coupling  coefficients. 

Considering  the  ka.,.!  eigenvalues  as  a  discrete  (infinite)  set  (k^i),  as  was  implied  for  ka 
by  the  use  ot  k^"),  then  we  can  define  the  discrete  coupling  coefficients  as 

ci"’  (a  +  1 ;  a) .  .jjL  [<(««)v,W(ki'ii) .  .  «-V,(«(ki»>)>„ 

(37) 

and  (36)  becomes: 

(ki"))  =  £  i  (a+1 ;  a)  (k^l^i)  ;  s  =  1, 2  (38) 

i  p»i 

The  coupling  coefficients  can  be  expressed  in  more  detail  when  the  specific  functional 
forms  of  the  eigenfunctions  appearing  in  the  inner  products  are  used  in  (37).  In  this  case  we  can 
use  the  orthogonality  of  the  vector  cylindencal  harmomes  to  reduce  the  coupling  factors  to  sim¬ 
ple  integrals  over  the  vertical  (z)  coordinate  on  the  boundanes  of  each  zone  Va.  Specifically, 
from  (37)  for  the  P-SV  case,  using  the  eigenfunction  expressions  given  earlier  in  (13)  -  (14),  one 
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has: 


RC,r^  (a+l;  a)  =  -^ 

+  |<(“+»)E, .  («)Sn>  -  ,  WEn>[B^f)(kA!ijp„)-B»^*)(ki'')pa) 


<(a+l)D, ,  (a)R^>  -  <(a+l)R,  ^  («)d^>  ^  >^)(k^'iiP«)-P^*)(k4")pa) 

(39) 


where  Oa  =  2npa,  with  po  denoting  the  value  of  the  radial  coordinate  on  the  surface  Ea-  Funher 
the  various  inner  products  involve  the  "vertical  eigenfunctions"  defined  in  (13)  :nd  (17);  where 
these  inner  products  have  explicit  forms  of  the  type: 


«"»D, .  i«R„>  ■  J  ;  ki''’W2. 


m) 


with  similar  expressions  for  the  other  products  in  (39).  If  these  products  are  compared  to  those 
in  (26)  and  (27)  •  or  more  directly  to  the  orthogonality  relations  involving  the  vertical  eigenfunc* 
tions  given  by  the  equations  following  equation  (27)  •  it  can  be  seen  that  the  inner  producu  in 
(39)  reduce  to  delta  functions  (^the  eigenfunctions  in  the  zones  Vq  and  Vq^.!  are  the  same:  that 
is.  if  s  ctc-  This,  of  course,  is  as  it  must  be.  since  only  when  the 

physical  prq>erties  in  the  two  zones  are  identical  will  ihe  eigenfunctions  be  the  same  and  it  then 
follows  that  the  coupling  matrix  must  be  diagonal  -  that  is  that  the  boundary  between  the  two 
zones  produces  no  vioss  mode  excitation  and  is  transparent.  We  see,  therefore,  that  t  le  analyti¬ 
cal  expression  in  (39)  for  the  coupling  does  indeed  have  this  required  property. 

The  normalization  factor  for  C  is  the  rano  rNs^®>  /  na  which  can  be  redefined  as  kNs^“’, 
where  from  the  previous  expressions  for  rNs^®),  in  (28)  and  (29),  this  constant  has  the  form: 


rN*(»>  = 


(40) 


Pa-l 


Pa 


(kit">  Pa-l)-PI^W9a-l)  +  B(?>(k4n>p„.,)-B*(f'(k4">Pa.l)H 
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In  an  exactly  analogous  fashion  the  coupling  coefficients  for  the  SH  modes  are  found  to  be; 


(a+1 a)  ^  [<(«+»F, .  «^)Tn>  - 

(41) 

where  the  inner  products  are  again  of  the  simple  form; 

<(a+i)F, .  (a)T„>  ,  ;  ki!i,)(«>T„(Zo  ;  ki"))  dZo  {41a) 

Further,  we  can  again  define  a  new  normalization  factor  /  n„  which  has  the  foi  m: 

•  /  \ 

lS'"*"  C'fl(ki">p.)CS;>(ki">P«)+  ^  Cfi>(ki">Pa-i)  Qif>(k4«P.-i)  (42) 

i  I  J  J 

The  corr  outarions  involved  in  determining  these  coefficients  are  straightforward,  since  the 
cylinderical  harmonics  are  tabulated  and  the  integrals  over  the  vertical  coordinate  Zo  involve 
simple  integn  Is  of  exponentials  that  can  be  evaluated  analytically,  in  closed  form,  for  the  gen¬ 
eral  case. 

Since  (3  0  constitutes  a  set  of  two  equations  for  s  s  l  and  s  s  2,  corresponding  to  forward 
and  backward  propaganng  modes  and  since  the  sums  on  the  right  can  clearly  be  expressed  as  a 
product  of  maihces,  it  is  natural  to  write  the  results  in  matrix  form.  Therefore  we  define; 

(k4‘>) 

m  ;  for  s  =  1  and  2  (43a) 

(k^N)) 

k  * 

and  a  similar  column  matnx  of  length  (L)  denoted  where  ihe  angular  index  m  has 

been  'luppressed  m  w’nnng  the  mode  excuation  matrices.  Further,  we  can  define  coupling 
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mrtrices  by: 


Cn<PC2i<P'*^  Cli^P'*) 
C,2(P**)  C22^P‘*’  Cl2<P'®^ 


CiN^P'*)  ...Cln^p-*) 


(43b) 


for  each  s  and  p  value,  where  s «  1,2  and  p  »  1,2.  With  these  definitions  one  c:  n  write  the  sys¬ 
tem  of  equations  implied  by  (38)  in  the  form: 


*  ' 

((«a«>l 

»  ' 

(WflOl) 

»  , 

.  ( 

[(a+t)flW] 

>  • 

where  the  forward  and  backward  propagating  mode  excitation  coefficients  are  •  tsplayed  expli¬ 
citly.  In  defining  the  matrices,  and  in  writing  the  matrix  result  m  (44),  the  "a  indices" 
have  been  suppressed.  However,  when  confusion  can  arise  they  should  be  written  as 
[C|n^P**>  (0  4- 1 :  a)],  etc.,  since  the  a  indices  change  when  the  matrix  refers  to  a  boundary  other 
than  £a'  (cf-  Between  the  zones  Va.i  and  V^,  on  the  surface  Ia-i>  the  cotpling  matrix  is 
expressed  as  (a ;  a- 1)1). 


Obviously  the  coupling  matrices  are  square  only  if  L  ■  N,  that  is  if  we  use  is  many  modes 
in  Vq  as  in  V«h,|  to  represent  the  propagating  waves.  This  choice  will  be  adhe:  jd  to,  from  this 
point  forward,  although  it  is  not  a  necessary  condmon. 

It  is  evident  that  the  partitioned  matrices  can  te  wntten  in  unpanittone  l  form  as  well, 
where,  with  L  >  N,  the  mode  coefficient  matrices  are  of  d  nension  (2.N  x  1)  and  the  coupling 
matrix  is  square  and  of  dimension  (2N  x  2N).  Thus,  we  can  also  define  mode  coefficient 
matrices  consisting  of  the  (ordered)  mode  coefficients  for  the  forward  and  backward  propagating 
modes  in  the  zones  Vq  and  Va^.}  as  (say): 
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[(a+Dfl/*)] 


and,  similarly,  we  can  define  what  can  appropnately  be  called  a  horizontal  propagator  matiix: 


A/in(a+l ;  a) 


1  L  J  L 


cm 


Now  the  equation  (44)  can  be  written  in  the  more  compact  form: 


(«>mn  *  //|n(a+l ;  a) 


and  expresses  the  required  conditions  between  the  mode  coefficients  in  neighboring  zones. 
If  we  take  successive  values  of  a,  with  a  ranging  from  1  to  M«1  say,  then  we  get 

[«>m„]«l/yu,(3;2))pfW,l 


Clearly,  by  noting  in  these  equations  that  the  indices  1  and  n  are  just  dummy  indices  providing  a 
numbering  system  for  the  eigenvalues,  then 

=  [//i„(2;l)]  [/yi„(3;2)]  ■■•[// i„(M;M-l)]  [w)/n, 

by  successive  subsmunons.  Consequently,  we  can  wnte,  for  any  p  ^  a  +  1: 
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This  is  a  propagator  equation  that  connects  the  mode  coefficients  in  any  zone  Vp  with  those  in 
any  other  zone  Vo.  In  case  P  =  a  +  1  the  equation  (47)  reduces  to  equation  (46),  which  connects 

the  coefficients  in  any  two  neighboring  zones.  Since  the  coupling  coefficients  composing  !  Hin  I 

L  j 

can  be  computed  from  the  simple  eigenfunction  inner  products  at  the  zone  interfaces,  this  equa¬ 
tion  provides  the  means  of  computing  mode  coefficients  that  produce  displacements  and  trac¬ 
tions  satisfying  all  the  boundary  conditions  along  the  venical  boundaries  of  the  medium.  Since 
the  eigenfunctions  used  already  satisfy  the  boundary  conditions  along  the  horizontal  boundanes 
in  each  zone,  then  by  use  of  the  horizontal  propagator  relation  all  the  boundary  conditions  in  tlie 
laterally  and  vertically  "layered"  medium  being  considered  can  be  sansfied. 

Summary  and  Conclusions 

The  basic  method  described  here  makes  use  of  normal  mode  expansions  of  the  wave  Held 
in  each  partitioned  sub-region  of  the  medium  within  which  the  medium  is  uniform  in  the  lateral 
directions.  Thus  the  medium  is  partitioned  into  laterally  uniform  zones  and  complete  normal 
mode  solutions  are  obtained  fro  each  horizontally  layered  zone.  In  the  analytical  development 
the  "zonal  eigenvalues  and  eigenfunctions"  are  generated  by  treating  each  zone  as  a  layered  half 
space  or  radially  layered  sphere,  as  ia  appropriate  for  th  -  medium  geometiy.  The  resulting  set  of 
modes  are  then  used  as  bases  for  expansions  of  the  w.tve  fields  in  the  layeied  subicgions.  Ilic 
mode  expansions  defined  on  the  zones  are  then  "connected"  by  matching  tequating)  the  exact 
Green’s  function  representations  of  the  wave  fields  in  each  zone  at  the  common  boundaries 
between  the  zones  where  contmuity  of  displacement  and  traction  is  required.  This  results  in  the 
defimtion  of  a  "lateral  propagator"  of  the  wave  field  when  applied  to  all  the  zones  mak  ng  up  the 
entire  medium  and  is,  in  application,  very  similar  to  the  classical  "venical  propagator  method. 
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The  method  is  exact  when  the  lateral  vanations  are  actually  discontinuous  step  changes  in  pro- 
pemes.  When  the  actual  changes  can  be  approximated  as  a  sequence  of  step  the  method  should 
be  superior  in  computational  accuracy  and  speed  to  numencal  methods. 

In  implementations  of  this  method  a  is  only  necessai7  to  compute  the  "zonal"  normal 
modes  once,  and  subsequently  these  zonal  mode  solutions  can  be  combined  in  a  variety  of  ways, 
using  the  lateral  propagator  equation,  to  produce  theoretically  predicted  wave  fields  in  many  dif¬ 
ferent  laterally  varying  structures  wintout  the  necessity  of  a  complete  recomputation  of  wave 
fields  in  each  new  structure.  Further,  the  propagators  are  analytically  defined  so  that  manipula¬ 
tions  related  to  inversion  and  perturbation  calculations  can  be  considered.  For  these  reasons,  and 
because  of  its  inherent  high  accuracy,  this  method  should  prove  useful  in  modeling  seismic  wave' 
fields  in  complex  media  and  in  inversion  studies.  In  the  present  study  the  method  is  developed 
in  detail  for  two  dimensionally  vanable  media,  using  cylindrical  coordinates  and  wave  functioas. 
However,  analogous  results  in  rectangular  and  spherical  coordinates  may  be  obtained  using  the 
same  procedure  and  are  app.''opnate  for  media  with  variability  in  all  three  spatial  dimensions. 
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Dr.  Dale  Glover 

DIA/DT-IB 

Washington,  DC  20301 

Mr.  James  C.  Batds 

Phillips  Laboratory/LWH 

HanscomAFB.1^  01731-5000 

Dr.  T.  Hanks 

uses 

Natl  Earthquake  Research  Center 
345  Middle&ld  Road 

Menlo  Park.  CA  94025 

Harley  Benz 

U.S.  Geolofteal  Survey,  MS-977 
34SMiddleliieldRd. 

Menlo  Park,  CA  9402S 

Dr.  Roger  Hansen 

AFTAC/rr 

Patrick  AFB,  FL  32925 

Dr.  Robert  Blartdford 

AFTAC/TT  • 

Center  for  Seismic  Studies 

1300  North  17th  Sl  Suite  1450 

Arlington,  VA  22209-2308 

Paul  Johnson 

ESS-4,  Mail  Stop  J979 

Los  Alamos  Natiooal  Laboratory  ' 
Los  Alamos,  NM  87545 

EricChael 

Division  9241 

SandiaUboitiocy 

Albuquerque,  NM  87185 

Janet  Johnston 

Phillips  Labonuory/LWH 
HanscomAFB,!^  01731-5000 

Dr.  John  J.  Cipar 

Phillips  Laboratorv/LWH 

Hans«:om  AFB,  \iA  01731-5000 

Dr.  Katharine  Kadinsky-Cade 
Phillips  Lab<mtoty/LV^ 
HanscomAFB.htA  01731-5000 

Cecil  Davis 

Group  P-15.  Mail  Stem  D406 

P.O.  Box  1663 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87544 

Ms.  Ann  Kerr 

IGPP,  A-025 

Setipps  Institute  of  Oceanogn^hy 
Univenity  oi  California,  San  Diego 
U  Jolla,  CA  92093 

Mr.  Jeff  Duncan 

Office  of  Congressman  Markey 

2133  Rayburn  House  Bldg. 

Washington,  DC  20515 

Dr.  Max  Kooniz 

US  Dept  of  Energy/DP  5 

Forres)^  Building 

1000  Independence  .Avenue 
Washington,  DC  20585 

Dr.  Jack  Evemden 
uses  •  Earthquake  Studies 

345  MiddlefieU  Road 

Menlo  Park.  CA  94025 

Dr.  Wil.K.  Lee 

Office  of  Earthquakes.  Volcanoes. 
AEngineeriog 

345  .Viiddlefield  Road 

Menlo  Park.  CA  94025 

A."!  Frar.ke! 

uses 

922  National  Center 

Reitcn,  VA  22G92 

Dr.  William  Leith 

U.S.  Geological  Survey 

Moil  Stop  928 

Reston,  Va  22092 
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Dr.  Richard  Lewis 

Director,  Earthquake  Engineering  &  Geophysics 

U.S.  Army  Corps  of  Engineers 

Box  631 

Vicksburg,  MS  39180 

Mr.  Chris  Paine 

Office  of  Senator  Kennedy 

SR  315 

United  States  Senate 
Washington,  DC  20510 

James  F.  Lewkowicz 

Phillips  Laboratory/LWH 

Hanscom.AFB,NlA  01731-5000 

Colonel  Jerry  J.  Perrizo 
AFOSR/NP,  Building  410 
Bolling  ATO 

Washington,  DC  20332-6448 

Mr.  AlftedUeberman 

ACDA/VI-OA'State  Department  Bldg 

Room  5726 

320 -2lst  Street,  NW 

Washington,  DC  20451 

Dr.  Prank  F.  Pilotte 
HQAFTAC/rT 

Patrick  AFB,FL  32925-6001 

Stephen  Mangbo 

Phillips  Laborittory/LWH 

HanscomAPB,MA  01731-5000 

Kade  Poley 

OA-AQSmdC 

Room  4X16NHB 

Washington,  DC  20505 

Dr.  Roben  Masse 

Box  25046,  Mail  Stop  967 

Denver  Federal  Center 

Denver.  CO  80225 

Mr.  Jack  Rachlin 

U.S.  Geological  Survey 
Geology,  Rm  3  C136 

Mail  Siw  928  Nadonal  Center 
Reston,VA  22092 

ArtMcGarr 

U.S.  Geological  Survey,  MS-977 

345  Middlefteld  Road 

Dr.  Robert  Reinke 

WUNTESG 

KirtlandAFB.NM  87117-6008 

Menlo  Park,  CA  94025 


Richard  Morrow 
AQ)A/VI,  Room  5741 
320  21st  Street  N.W 
Washington,  DC  2045 1 


Dr.  Carl  Newton 

Los  Alamos  National  Laboratory 

P.O.  3ox  1663 

Mail  Stop  C335.  Group  ESS-3 
Los  Alamos,  NM  87545 

Dr.  l*ao  Nguyen 
AFTaQTTR 
Patrick. AFB.FL  32925 


Dr.  KiP.r.eth  H.  Olsen 

Los  .ALunos  Sciendrlc  Labontorv 

P.  O.  Bo.x  1663 

Mail  Stop  D-iC6 

Los  Alamos,  NM  87545 


Dr.  Byron  Ristvet 

HQ  DNA,  Nevada  Operadons  Office 
Atm:  NVCG 
P.O.  Box  98539 
Las  Vegas,  NV  89193 

Dr.  George  Rothe 
HQAFTAC/rrR 
Patrick  AFB,FL  32925-6001 


Dr.  Alan  S.  Ryall,  Jr. 
DARPA/NMRO 
1400  Wilson  Boulevard 
ArUngton,  VA  22209-2308 


Dr.  Michael  Shore 
Defease  Nuclear  .Agency/SPSS 
6801  Telegraph  Road 
.Mexandria,  W.\  22310 


Phillips  Laboratory 
Atm:  XO 

Hanscom  AFB,  MA  01731-5000 


Dr.  Laity  Turnbull 
ClA-OSWR/J^ 
Washihgton,  DC  2050S 


Dr.  Thomas  Weaver 
Los  Alamos  National  Laborauxy 
F.O.  Box  1663,  MaU  Stop  C33S 
Los  Alamos.  NM  87S4S 


Phillips  Laboratory 
Research  Library 
ATTN:  SULL 

Hanscom  AFB .  MA  01731-3000 


Phillips  Laboratory 
ATTN*  SUL 

KiitlaiMlAFB.NM  87117-6008 


Secretary  of  the  Air  Force 
(SAFRO) 

Washington.  DC  20330 


OOice  of  the  Secretary  Defense 
DDR&E 

Washington.  DC  20330 


HQDNA 

Atm:  Technical  Library 
Washington.  DC  2030S 


DARPA/RMO^RETRIEVAL 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


DARP.X/'RMO/Secariiy  Office 
1400  Wilson  Boulevard 
.Arlington,  V A  22209 


l^ips  Laboratory 
Atm:  LW 

Hanscom  AFB.  MA  01731-5000 


DARPA/PM 

1400  Wilson  Boulevard 

Arlinron,VA  22209 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria.  VA  22314  (2  copies) 


Defense  Intelligence  Anncy 
Directomte  forSdendnc  &  Technical  Intelligence 
(2  eopits)  Attm  DTIB 

Washington.  DC  20340-6138 


AFTAOCA 

(STINFO) 

Patrick  AFB.  a  32923-6001 


TACTEC 

Batielle  Memorial  Institute 
305  King  Avenue 

Columbus,  OH  43201  (Final  Repon  Only) 
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CONTRACTORS  (FOREIGN) 


Dr.  Ramon  Cabre,  S.J. 
Observatorio  San  (jalixto 
asiUa5939 
La  Paz,  Bolivia 


Prof.  Hans-Petcr  Haijes 
"  Institute  for  Gcophysik 
Ruhr  Univcrsity/Bochum 
P.O.  Box  102148 
'  4630  Bochum  l.FRG 

Prof.  Eystein  Husebye 
NTNF/NORSAR 
P.O.  Box  SI 

N-2007KjcUcr.  NORWAY 


Prof.  Brian  LN.  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
GP.O.  Box  4 

Canberra  2601,  AUSTRALIA 

Dr.  Bernard  Massinon 

Soclete  Radiomana 

27  rue  Claude  Bernard 

75005  Paris,  FRANCE  (2  Copies) 


Dr.  Pierre  .Mechelcr 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris,  FRANCE 


Dr.  Svein  Mykkeitveit 
NTNF/NORSAR 
P.O.  Box  51 

N-2007Kjeller,NORW.AY  (3  copies) 
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FOREIGH 

Dr.  Peter  Basham 

Earth  Physics  Branch 

Geologit^  Survey  of  Canada 

1  Observatory  Crescent 

Ottawa.  Ontario,  CANADA  K1 A  0Y3 

Dr.  Eduard  Berg 
Institute  of  Geophysics 
University  of  Hawaii 
Honolulu.  HI  96822 


Dr.  Michel  Bouchon 
I.R,I.G.M.-B.P.  68 
38402  St.  Martin  D'Heres 
Cedex,  FRANCE 


Dr.HilmarBupgum 
NTNF/NORSAR 
P.O.  Box  51 

N-2(X)7KjcUer,  NORWAY 


Dr.  Michel  CampUlo 
Observatoire  de  Grenoble 
I.R.I.G.M.-B.P.  53 
38041  Grenoble,  FRANCE 


Dr.  Kin  Yip  Chun 
Geophysics  Division 
Physics  Department 
University  of  Toronto 
Ontario.  CANADA  M5S  1A7 

Dr.  Alan  Douglas 
Ministry  of  Defense 
Blacknest,  Brimpton 

Reading  RG7-4RS,  UNITED  KINGDOM 


(OTHERS) 

Dr.  Tormod  Kvaema 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Peter  Marshall 
Procurement  Executive 
Ministty  of  Defense 
Blacknest.  Brimpton 

Reading  FG7-4RS.  UNITED  KINGDO.VI 

Pre  f.  Ari  Ben-Menahem 
Department  of  Applied  Mathemadcs 
Weizman  Institute  of  Science 
Rehovot,  ISRAEL  951729 


Dr.  Robert  North 

Geophysics  Division 

Geolo^cal  Survey  of  (Canada 

1  Obs^atoty  Cr^cent 

Ottawa,  Ontario,  CANADA  KIA  0Y3 

Dr.  Erode  Ringdal 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller.  NORWAY 


Dr.  Jorg  Schlittenbardt 

Federal  Institute  for  Geosciences  &  NaiT  Res. 
Postfach  510153 

D-3000  Hannover  51,  FEDER.AL  REPUBL  iC  OF 
(jERMANY 


Universita  Degli  Studi  Di  Trieste 
Facolta  Di  Ingegneria 

Istituto  Di  Mmiere  E.  Geofisica  Applicata,  Trieste, 
ITALY 


Dr.  Manfred  Hengcr 

Federal  Institute  for  Geosciences  &  Nat!  Res. 

Postfach  510153  Dr.  John  Woodhoi  se 

D-3000  Hanover  51.  FRG  Oxford  University 

Dept  of  Earth  Sciences 
Parks  Road 

Ms.  Eva  Johannisson  Oxford  0X13PR,  ENGLAND 

Senior  Research  Officer 
National  Defense  Research  Inst. 

P.O.  Box  27322 
S-102  54  Stockholm,  SWEDEN 

Dr.  Fekadu  Kebede 

Geophysical  Observatoiy,  Science  Faculty 
Addis  .Ababa  University 
P.  0.  Bo.x  1176 
Addis  Ababa,  ETHIOPIA 
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